Nonlinearities in modified gravity cosmology — II. Impacts of modified gravity on the 

halo properties 
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The statistics of dark matter halos is an essential component of understanding the nonlinear 
evolution in modified gravity cosmology. Based on a series of modified gravity N-body simulations, 
we investigate the halo mass function, concentration and bias. We model the impact of modified 
gravity by a single parameter f, which determines the enhancement of particle acceleration with 
respect to GR, given the identical mass distribution (£ = f in GR). We select snapshot redshifts such 
that the linear matter power spectra of different gravity models are identical, in order to isolate the 
impact of gravity beyond modifying the linear growth rate. At the baseline redshift corresponding 
to zs = 1.2 in the standard ACDM, for a 10% deviation from GR(|£ — 1| = 0.1), the measured halo 
mass function can differ by about 5 — 10%, the halo concentration by about 10 — 20%, while the 
halo bias differs significantly less. These results demonstrate that the halo mass function and/or 
the halo concentration are sensitive to the nature of gravity and may be used to make interesting 
constraints along this line. 
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PACS numbers: 98.80.-k; 98.65.Dx; 95.36. 



I. INTRODUCTION 

Dark matter halos are prominent structures in the dark 
universe. Their abundance, density profile and clustering 
(halo bias) contain valuable cosmological information. In 
particular, dark matter halos form and grow under grav- 
itational instability. Hence halo properties contain rich 
information on the nature of gravity at ~ Mpc scales 
and above, and can provide strong gravity constraints 
These properties also provide a powerful tool to 
understand the matter clustering through the halo model 
[IH [lfjl], which is also a sensitive measure of gravity. 

In [20] (hereafter paper I) we run a controlled set of N- 
body simulations with identical initial condition to study 
the impact of modified gravity (MG) on the matter power 
spectrum. These simulations adopt a MG parameteriza- 
tion with a single parameter £. £ is the relative enhance- 
ment of nonrelativistic particle acceleration with respect 
to general relativity (GR), given the identical mass dis- 
tribution. More specifically, it is the quantity that enters 
into the vb-p relation in Fourier space, 



k 2 yb = -(4:irGa 2 6p 



(1) 



Here, dp = p — p and p is the matter density. i\) is defined 
through ds 2 = -(l+2y!j)dt 2 +a 2 (l + 2(f>)dx 2 . vb is the only 
gravitational potential that nonrelativistic particles can 
sense. So the value of £(k, z) as a function of scale and 
redshift, along with the initial condition and the expan- 
sion rate of the universe, completely fixes the evolution of 
the matter density and velocity fields. This is the major 
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reason that we adopt this single parameter parametriza- 
tion on modified gravity N-body simulations. For more 
detailed discussion on this parametrization, please refer 
to paper I. 

Neverthless, such parameterization is highly simplified 
in the sense that it lacks any screening mechanism [l?} 
required to pass the solar system tests. For this reason, 
these simulations lose the capability to explore the rich 
consequences induced by the MG environmental depen- 
dence. These features have been explored in advanced 
simulations on f(R) and DGP [H, f2M2^ . Nevertheless, 
our simulations benefit from costing no extra time than 
the ordinary CDM simulations. So in principle one can 
run a large number of such simulations within reason- 
able amount of time to fairly sample relevant parameter 
space. One then interpolates/extrapolates the simulation 
results to explore the whole relevant parameter space and 
understand the nonlinear evolution for £ of general spa- 
tial and time dependence. The hope is that, by choose 
appropriate spatial and time dependence in £, one can ef- 
fectively take the environmental dependence into account 
& 

In the current paper, we will examine the halo prop- 
erties in these simplified MG simulations, including the 
halo abundance, concentration and bias as a function of 
mass and redshift. In the standard framework of struc- 
ture formation in GR, these properties, in particular the 
halo abundance and bias, are largely fixed by the linear 
matter power spectrum at the corresponding epoch. We 
utilize this property to better isolate the impact of MG. 
We will compare the above properties not at the same 
redshift, but at redshifts of identical linear matter power 
spectrum. The same method is also adopted in paper I. 
It has extra benefit of reducing cosmic variance and nu- 
merical artifact, since we focus on the differences between 
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different MG simulations with identical initial condition. 

This paper is organized as follows. We briefly describe 
our simulations in $Il]and present results in §1111 We dis- 
cuss and summarize in §JV] 

II. N-BODY SIMULATION 

We run a set of simulations with the TreePM paral- 
lel code GADGET-2 H3] at Shanghai Astronomical Ob- 
servatory. All the simulations evolved 512 3 dark mat- 
ter particles in a periodic box of 300/i _1 Mpc on a side. 
The cosmological parameters used in the simulations are 
= Hdm + = 0.276, fl h = 0.046, h = 0.703, 
tt A = 0.724, n = 1, and cr 8 = 0.811. The particle 
mass and softening length are 1.541 x 10 10 /i _1 M and 
12.89 kpc, respectively. Glass-like cosmological initial 
conditions were generated at redshift z — 100 using the 
Zel'dovich approximation. 

In this paper as in Paper I, the modified gravity model 
is characterized by a single parameter £, which deter- 
mines the enhancement of particle acceleration with re- 
spect to GR, given the identical mass distribution (£ = 1 
in GR). All the simulations for different £ values have the 
same linear power spectrum, such that 

P L (k;z Sl (=l)=P L (k;z (), (2) 

where Pl is linear matter power spectrum, zs and z^ 
denote the redshifts in the standard ACDM and the MG 
cosmology, respectively. Since all the simulations begin 
with the identical initial condition, the above relation 
means that 

D(z s ,C = l)=D(z c ,Q, (3) 

where D(z, £) is the linear density growth factor. Given 
redshifts z$ in the ACDM, through Eq. ^ we can cal- 
culate the corresponding redshifts which are shown 
in Table I of Paper I. In this paper, we mainly focus on 
the redshift z$ = 1.2 in the ACDM and the correspond- 
ing redshifts in the MG universe (10th row of Table I in 
Paper I). 

Dark matter halos were identified from the simulation 
at each snapshot using the standard friends-of-friends 
(FOF) algorithm^ with a linking length of b = 0.2 
times the mean interparticle separation. For currently 
favored cosmologies, 15-20% of b = 0.2 FOF halos have 
irregular substructure or have two or more major halo 
components linked together 29] . Only isolated, relaxed 
halos are well-fit by the Navarro-Frenk- White (NFW) 
profile [30|, therefore, we use the virial halo mass to fit 
NFW profile in the simulations. The mass M200 is 
defined as the interior mass within a sphere of radius 
-R200 centered on the most bound particle of the halo, 
where -R200 is the radius within which the over-density is 
A200 = 200 times the mean density /3 moa n of the universe. 
Thus, the mass and radius of a halo are related by 

M200 = o^^aOOPmcan^OO' ( 4 ) 
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FIG. 1: Comparison of halo distributions from N-body simu- 
lations for different £ values at the baseline redshift zs = 1.2, 
which is that of £ = 1 (ACDM). The panels show the dark 
matter halos with the mass M 200 > 10 11 ' 5 /i _1 M Q in slices 
100 x 100 x 8/i _1 Mpc, and the sizes of the dots are propor- 
tional to the radii -R200 of the halos. 



where p me an is the mean density of the universe at red- 
shift z. Note that some authors use a value A v j r moti- 
vated by the spherical collapse model, where the corre- 
sponding radius is specified as i? v i r . This value is both 
redshift and gravity dependent. It gives 

^^vir — 339.5 

at z = and 178 at z ^> 1 for our adopted ACDM pa- 
rameters |3ll|. It also shows significant deviation in MG 
models [14j |. Since we will compare the halo properties 
of different MG models at different redshifts (e.g. z s vs. 
Z(), this latter choice could mix the impact of background 
cosmology with the impact of gravity. For this reason, 
we will not adopt this latter definition. Unless other- 
wise specified, throughout the paper we use a fixed value 
A = 200 to define the halo mass. 

Fig. Q] shows the distributions of dark matter halos 
with the mass M200 > 10 115 /i _1 M Q in a slice of thick- 
ness 8ft. _1 Mpc for different £ values at the baseline red- 
shift z s = 1.2, which is that of ( = 1 (ACDM). Accord- 
ing to visual inspection, the general appearances of the 
large-scale structures are remarkably similar, as a result 
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FIG. 2: Measured halo mass functions from N-body simu- 
lations for different £ values at redshift zq where the linear 
matter power spectrum is identical to that in ACDM (£ = 1) 
at the baseline redshift zs = 1.2. Since all simulations begin 
with the identical condition, the measured halo mass function 
has similar Poisson error. Hence we only show the Poisson 
error bars of ACDM to demonstrate the measurement uncer- 
tainty. The lower panel shows the ratio of the measured mass 
functions of MG models to the one in ACDM when they have 
the same linear power spectrum. The plotted errors are the 
same Poisson errors showing in the upper panel. Since the 
Poisson errors of different simulations have strong positive 
correlations, they largely cancel when taking the ratio. So 
the plotted errors represent a conservative upper limit of the 
halo mass function ratio. The lower panel then shows that 
the observed deviations from a universal mass function are 
significant. 



of the same linear power spectrum for different £ val- 
ues. Although the shape of Cosmic Web (the pattern on 
clusters, filaments, sheets, and voids) is more or less pre- 
served in the £ ^ 1 runs when compared to the ACDM 
case, we can depict arising differences when looking into 
the internal properties of dark matter halos. 



III. NUMERICAL RESULTS 



In this section we show the halo mass function, concen- 
tration and bias as a function of M2oo- Structure grows 
more slowly in £ < 1 models than that in GR. So for too 
low zs in ACDM, structure growth in a £ < 1 universe 
may not be able to catch up with that of ACDM even at 

= (table I, paper I). For this reason, we only com- 
pare the results corresponding to Zg = 1.2 in ACDM, for 
which we can compare the halo statistics for all £. 



A. The halo mass function 

One of the long-standing efforts in precision cosmology 
is to determine the mass function of dark matter halos 
dn/dM, which is the number of halos per unit volume 
per unit mass. It is a sensitive probe of gravity and has 
been used to pu t strong constraints on modified gravity 
models (e.g. [24( on f(R) gravity). 

Although the number density of halos of a given mass 
depends on the shape and amplitude of the power spec- 
trum, analytical work has suggested that the halo abun- 
dance can be expressed by a universal functional form 
when expressed in terms of suitable variables [33 - l33 | . A 
convenient form to describe halo abundance can be ex- 
pressed as [35[ 



dM y ' ' 



n a >T7 777 » (°) 



M 



dM 



where po is the background matter density at redshift 
z = 0, and u 2 (M, z) is the variance of the linear matter 
power spectrum over a length R, 



a 2 (M,z) 



2tt 2 



k 2 P(k)W 2 (kR(M))dk, (6) 



where W(x) = 3[sin(x)— xcos(a;)]/a; 3 is the Fourier trans- 
formation of the top-hat filter, R(M) = (3M/4-Kp ) 1/3 is 
the smoothed radius with halo mass M, and D(z) is the 
growth factor. 

This definition of mass function has been widely exam- 
ined against N-body simulations and useful fitting func- 
tions f(a) have been provided by several authors [35l438| . 
Nevertheless, deviations from the universality have been 
reported. For example, Bhattacharya et al. [38] found 
through their N-body simulations that, for their wCDM 
cosmological models (where the dark energy equation of 
state parameter u> is constant in time, but uj ^ —1), the 
universality of the mass function is systematically broken 
at a level of 10%. Deviations from the universality are 
also detected in interacting dark energy models Cui et al. 
[39[ | . Cui et al. [39| reported that deviations can exceed 
~ 10% for most of the models in the high mass end. 

Fig. [2] shows the measured halo mass functions from 
our simulations for different £ values at the baseline red- 
shift z s = 1.2, which is that of £ = 1 (ACDM). The 
M 2 oo halo mass is used in the plot in order to get con- 
sistent comparisons with concentration and bias, where 
the same mass definition are adopted. We have checked 
our measured halo mass functions for ACDM simulation, 
which agree well with the fitting function proposed by 

m 

In order to assess the difference, the Poisson error bars 
are only added for £ = 1. Because of the same initial 
condition been used, the error bars for £ ^ 1 have the 
similar forms. In Fig. [2] the lower panel shows the ratio of 
the measured mass functions of MG models to the one in 
ACDM. The Poisson error bars for £ = 1 are added to as- 
sess the maximum limit. Since D(zq,Q = D(z$,C = 1); 
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the linear power spectrum and hence cr(M, z) are iden- 
tical. If the halo mass function is indeed universal and 
described by Eq. (O and Eq. (j6|), the mass function in 
the MG model should be the same as that in ACDM. 
However, the result in Fig. [2] shows the difference of 
mass functions between £ = 0.8 and C = l.O(ACDM) is 
about 20% with the halo mass range from 10 115 ft _1 M Q 
to lQ xw h~ x MQ. The difference between C = 1-2 and 
C = 1.0 is about 10%. We find that the difference of 
( = 1.2 is about twice that of ( = 0.8, although they 
have the same 20% deviation from GR. For a 10% devi- 
ation from GR(£ = 0.9, or £ = 1.1), the measured mass 
functions can differ by about 5 — 10%. Interestingly, the 
difference of C = 0.9 is also about twice that of ( = 1.1. 
In general, for the identical deviation of £ from unity, the 
deviation of halo mass function of £ < 1 is about twice 
that of C > 1, which implies that the decrease of gravity 
can lead to more difference of halo mass function. This 
behavior can also be found in the bottom-right panel of 
Fig. 4 in Paper I for the comparison of the nonlinear 
power spectra. For the same deviation of £ from unity, 
the deviation of nonlinear power spectra of £ < 1 are 
about twice that of £ > 1 at A 2 ~ 10. 

If the halo abundance is completely determined by the 
shape and amplitude of the linear power spectrum, the 
halo mass functions for different £ simulation should have 
the identical form. However, based on the measured halo 
mass functions from N-body simulations for different £ 
values at the same linear power spectrum, we find that 
this is not the case, which means that the halo abun- 
dance also depends on the structure growth history and 
the underlying gravity. The deviation of halo abundance 
becomes larger if the deviation of £ from unity is larger. 
The halo mass function is therefore a sensitive probe for 
our modified gravity model. 



B. Halo concentration 

The concentration parameter has been widely used to 
describe the internal structure of halos [42lT4q . Recently, 
Jeeson-Daniel et al. [13] investigated the correlation be- 
tween nine different dark matter halo properties, and 
they claimed that while the scale of a halo is set by its 
mass, the concentration is the most fundamental prop- 
erty. In this paper, the NFW form is used to fit the halo 
density profiles, which can be approximated by a simple 
formula with two free parameters, 



P( r ) 

Pmcai 



(r/r s )(l+r/r s )^ 



(7) 



where r s is a scale radius and S c is a characteristic den- 
sity contrast. The concentration of a halo is defined as 
c = -R200 / r s , thus S c and concentration are linked by the 
relation 
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FIG. 3: Measured halo concentrations from N-body simula- 
tions for different £ values at the baseline redshift zs = 1.2. 
The lower panel shows the ratio of the measured halo concen- 
trations of MG models to the one in ACDM. For clarity, we 
plot the error bars of the £ — 1 result. As explained in Fig. 
El these errors are conservative upper limit. Hence the ob- 
served differences in c are statistically significant. Difference 
in the nature of gravity contributes to these differences. How- 
ever, since zs 7^ z^, difference in c also arises from difference 
in the background matter density. More systematical studies 
are required to disentangle the two and to infer the nature of 
gravity from the halo concentration. 



In our calculation of concentration from the simulation, 
the radius of a halo is divided in uniform logarithmic bins 
(Alog 10 = 0.1), starting from the radius where the bin 
contains at least 20 particles to i? max = (or to i? v i r 
in case halo masses are defined using the over-density 
A v ; r ). Due to the relation of equation [8l at given halo 
mass M200, there is a single free parameter in equation [3 
which can be expressed as the concentration parameter 
c. The best-fit concentration parameter can be computed 
from simulation by minimizing the rms deviation a be- 
tween the binned p(r) and the NFW profile, 



a' 



N bh 



E 



[logio Pi - log 10 Pnfw(c)] 



(9) 



200 



3 [ln(l + c) - c(l + c)] ' 



(8) 



Fig. [3] shows the measured concentration of the halos 
contained at least 200 particles for different £ values. In 
each mass bin, the halos are randomly divided into five 
parts. Then the average concentration of each part is 
calculated. The error bars plotted in Fig. [3] is the stan- 
dard deviation of five average concentrations. In general 
the errors of the average concentrations of the total ha- 
los can be reduced by a factor of \/5. For £ = 1, we find 
that the concentrations of halos with iV > 1000 parti- 
cles (M 2 oo > 13.18 h~ 1 M Q ) can be approximately fitted 
by a power law c oc M^j 1 , which is in good agreement 
with the results of Neto et al. [451. For the halos with 
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N < 1000 particles, we find the power law index is less 
than —0.1, which is in good agreement with the results 
calculated by the model of Zhao et al. [52^ . 

Because the main focus in this paper is to compare 
the difference of concentration for different £ values, in 
the lower panel of Fig. [3j the ratio of the measured halo 
concentrations is plotted. For a 10% deviation from GR 
(|C — 1| < 0.1), the deviation of halo concentration can 
differ by 10 — 20%, which means that the modified gravity 
can significantly affect the structure formation on small 
scale and the internal properties of the dark matter halos. 
For £ = 0.8, the deviation of halo concentration even 
differ by larger than 50% compared to that of £ = 1.0. 
Here again, with similar difference in £, the concentration 
difference of £ < 1 is much larger than that of £ > 1 . 

The discrepancy between the different £ values is 
mainly due to the different background density although 
they have the same linear matter perturbation. Since 
structure grows faster in £ > 1 cosmology {zq > zs), 
halos in this universe form in a background with higher 
mean density (oc (1 + z) 3 ). We then expect them to have 
a smaller concentration [46j. For the same reason, we ex- 
pect halos with £ < 1 have a larger concentration. The 
behavior of the concentrations for different £ values is 
roughly similar to the trend of the concentrations cor- 
responding to the redshift z^ in ACDM models. Large 
values of the ratio of the concentration parameters from 
unity implies that the halo concentration is a valuable 
property to detect different modified gravity models. 



C. Halo bias 

We define the bias of dark matter halos as the ratio of 
the halo mass cross-power spectrum to the dark matter 
power spectrum 

6(fc,M 200 ) = — p ttt — > (10) 

* ram y^) 

where P nm (k, M) denotes the cross-power spectrum with 
halos of mass M200, and P m m(k) is the dark matter power 
spectrum. This measure dose not require a shot-noise 
correction, and it yields better statistics when the halos 
become sparse. The power spectrum is calculated us- 
ing Daubechies wavelet mass assignment, which avoids 
the sampling effect. Although halo bias is scale depen- 
dent in the quasi-linear and nonlinear regime, here we 
focus on the large-scale bias, where b is independent 
of k. We calculate b as the average over the 5 largest 
wavelength modes (k < 0.1/i/Mpc) in the simulation. 
We also check these results against bias as defined by 
bhh = v Phh/Pmm, and we find that there is a good 
agreement between bhm and bhh- 

In the ACDM cosmological model, the linear halo bias 
can be expressed as a function of v, [41], l48Tj50| . where 
v = 5 c /a is the ratio of the critical over-density required 
for collapse to the rms density fluctuation. 
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FIG. 4: Measured halo bias from N-body simulations for dif- 
ferent £ values at the baseline redshift zs = 1.2. The lower 
panel shows the ratio of the measured halo bias of MG mod- 
els to the one in ACDM. As in previous plots, the errors are 
that of ACDM. As explained in previous plots, these errors 
are conservative upper limit. For this reason, the observed 
difference in the halo bias is statistically significant, despite 
the apparently large error bars. 

Fig. |4] shows measured halo biases from N-body sim- 
ulations for different £ values at the baseline redshift 
zs = 1.2. The error bars is the standard deviation of the 
5 largest wavelength modes. The error bars for different 
£ values have very similar amplitude and only £ = 1.0 
has been plotted in the figure. The lower panel shows 
the ratio of the measured halo bias of MG models to the 
one in ACDM. We find that the difference of halo bias for 
different C values becomes larger as the increase of halo 
mass. However, the deviations of halo bias differ by less 
than 10% from GR (£ — 1) for all the £ 7^ 1 values used in 
our simulation. This behavior is expected due to the halo 
bias is a function of v in ACDM model. The simulations 
for different £ values have the same linear power spec- 
trum and growth factor, therefore, according to Eq. ^ 
they have the same density fluctuation a. Besides S c is 
weakly dependent on the redshift, thus the difference of 
the halo bias for different £ values is very small. We can 
conclude that the halo bias is a weak statistics to detect 
the modified gravity model. 



IV. DISCUSSION AND CONCLUSION 

In this paper we have compared the halo mass func- 
tion, concentration and bias based on a series of modified 
gravity N-body simulations. The modified gravity model 
is characterized by a single parameter £, which deter- 
mines the enhancement of particle acceleration with re- 
spect to GR. All the simulations for different £ values are 
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FIG. 5: Same as Fig. [2] but for the virial halo mass M v i r 
defined by a value A v i r from the spherical collapse model. 



FIG. 7: Same as Fig. [4] but for the virial halo mass M v i r 
defined by a value A v i r from the spherical collapse model. 
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FIG. 6: Same as Fig. [3] but for the virial halo mass M v i r 
defined by a value A v i r from the spherical collapse model. 

started from the same initial condition. The rcdshifts for 
different £ of the comparison of halo properties are se- 
lected when they have the same linear power spectrum 
as ACDM simulations at z$ = 1.2. 

Based on our various comparisons, we summarize our 
findings as follows. 

• The measured halo mass functions can differ from 
that in GR by about 5—10% for 10% deviation from 
GR (|C - 1| = 0.1). Although the difference must 



be symmetric with respect to £ = 1 in the limit 
( — » 1, it already shows significant asymmetry for 
|£— 1| ~ 10%. For example, the difference between 
the halo mass function of £ = 0.9 and that of GR 
is about twice of the corresponding case of £ = 1.1. 

• For a 10% deviation from GR(|C~1| = 0.1), the de- 
viation of halo concentration can differ by 10 — 20%, 
which shows that modified gravity can significantly 
affect the structure formation on small scale and 
the internal properties of the dark matter halos. 

• The halo bias is less sensitive to the nature of grav- 
ity. For all the C 7^ 1 values used in our simulation 
(0.8 < C < 1-5), the halo bias differ from that of 
GR by less than 10%. 

We note that all the above findings are obtained based 
on the halos that are defined to be the spherical over- 
density regions within which the average over-density are 
A200 = 200 times the mean density p mcan of the universe. 
This definition is widely used in literature and is easy to 
be implemented both in simulations and in observations. 
The impact of using different A has been discussed using 
high resolution N-body simulations [37j. 

One might argue that the differences between MG 
models and the GR shown above are mainly due to the 
definition of the halos. To test to what extent halo defini- 
tion impacts our results, we have also compared the halo 
mass function, concentration and bias using the virial 
halo mass M v i r (defined by a value A v j r [31| from the 
spherical model) instead of M200- According to the virial 
mass definition from GR, we do see in Fig. [5] that for a 
20% deviation from (|£ — 1| < 0.2), the deviation of halo 
mass function differs by less than 10%, which is smaller 
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than that of the halo mass M2oo- However, we can still 
detect the observed differences of the halo mass functions 
for different £ values. 

In Fig. O for ( > 1.0, the difference of c V i r (C) with 
respect to that of GR is roughly the same as the case 
of C2oo(C)- Since at redshift zs = 1.2, the universe is 
dominated by dark matter and i?200 is approximately 
equal to i? v ir, which is 195.7 for our adopted parameters. 
In one word, by comparing the results shown in Fig. [3] 
and Fig. [SJ we can see that there exists a significant dif- 
ference of the halo concentration for different C values, 
which illustrates that deviations from general relativity 
can strongly affect the mass distribution in the nonlinear 
regime, although the large-scale structures are remark- 
ably similar (see Fig.[I| [53j . 

Fig. [7] shows that the bias of halos defined with virial 
mass still shows weak dependence on gravity, for all the 
( ^ 1 considered (0.8 < ( < 1.5). This confirmes pre- 
vious finding that the halo bias may not be a sensitive 
probe of gravity. This conclusion would be robust for 
halos of mass less than 10 13 /i~ 1 M Q . Neverthless, our 
simulations do not have sufficiently large box size to ro- 
bustly measure the bias of halos with mass larger than 
1O 13 /i _1 M and evaluate its dependence on gravity. 

In th above results A v j r is evaluated assu ming GR. It 
is known that A v j r is changed in MG models |l4t l40j . To 
what extent it can account for the found difference in the 
mass function (l^.l40j is an issue for further investigation. 



Our results have robustly demonstrated that the halo 
mass function is not completely determined by the shape 
and amplitude of the linear power spectrum. As a con- 
sequence, it shows significant dependence on the nature 
of gravity. We also show that the halo concentration is 
also affected by gravity. These results can be used to 
understand the impact of gravity on halo formation and 
distribution. They can also be used in halo model to 
understand the nonlinear matter clustering (20j . 
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